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The central region of the dark matter halo of the Milky Way is a promising target for a search for a particle dark 
matter annihilation signal. The H.E.S.S. Collaboration has published a search for a photon flux originating from 
dark matter particles annihilating in the galactic center region. No significant excess was observed and upper 
limits on the velocity averaged dark matter self annihilation cross section were derived. The limits exclude 
the self annihilation of dark matter particles in the ~ 1 TeV to ~ 4 TeV mass range with a velocity averaged 
annihilation cross section larger than ~ 3 ■ 10~^^ cm'^/s, i.e. about one order of magnitude above the prediction 
for a thermal relic dark matter particle. A detailed and realistic Monte Carlo study of new strategies for the 
search for a particle dark matter annihilation signal from the Milky Way dark matter halo with the High Energy 
Stereoscopic System is presented and the sensitivity of different experimental approaches is compared. 
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1. Particle Dark Matter and the Milky Way 
Dark Matter Halo 

The dynamics within the Milky Way points towards 
the existence of a large amount of invisible, i.e. dark, 
matter in the Milky Way. Many extensions of the 
standard model of particle physics naturally predict 
the existence of yet unobserved weakly interacting 
massive particles (WIMPs) that do not directly couple 
to light but only to the weak sector of the standard 
model. It is thus inviting to explain the failure to un- 
derstand the dynamics within the Milky Way based 
on the visible distribution of matter by the additional 
abundance of WIMPs. The existence of WIMPs is also 
in agreement with cosmological structure formation 
arguments and precision measurements of the dark 
matter density in the universe. 

In many models, WIMPs are expected to be stable 
on cosmological timescales but have the ability to an- 
nihilate into standard model particles. In the mini- 
mal supersymmetric extension of the standard model 
of particle physics with conserved R-parity, the neu- 
tralino, i.e. the superposition of supersymmetric part- 
ners to the neutral electroweak gauge bosons, is a 
Majorana particle dark matter candidate and is thus 
even able to self annihilate. Depending on the an- 
nihilation channel, the standard model annihilation 
products lead in part to a continuous photon spec- 
trum. The photon energies can, depending on the 
particle dark matter mass, be in the Very High En- 
ergy (VHE, E > 100 GeV) 7-ray regime. The number 
of expected particle dark matter annihilation events 
in a given volume is proportional to the dark matter 
density squared. Large dark matter densities are thus 
promising for a search for dark matter self annihilation 
events. Recent large scale N-body simulations of grav- 
itationally interacting pseudo particles predict that 
the dark matter density distribution at the spatial 
scale of galaxies can be parametrized by simple func- 
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Figure 1: Two parametrizations of the dark matter 
density distribution as a function of the distance to the 
galactic center resulting from large scale N-body 
computer simulations. The two distributions correspond 
to the best fit to the dark matter distribution in two 
different Milky Way mass scale N-body comp uter 
simulations as obtained in iPieri et al. I 120 llH . The dark 
matter density distributions agree for distances larger 
than ~ 45 pc within a factor of two. Both distributions 
are normalized to a dark matter density of 0.3 GeV/cm'^ 
at a distance of 8.5 kpc to the galactic center, i.e. the 
radial distance of the sun to the galactic center. 



tions and few parameters (see f.i. iPieri et al. I (201lj ). 
Figure [T] shows the dark matter density distributions 
as recently obtained in two different Milk y Way mass 
scale N-body computer simulations (see Pieri et al. I 
j2011 |). Both distributions are normalized such that 
the density at the position of the sun agrees with the 
value of 0.3 GeV/cm'^ derived from different dynam- 
ical constraints. A central result of N-body simula- 
tions involving only dark matter particles is that the 
dark matter density is in general steeply increasing to- 
wards the center of galaxies. The agreement between 
the resulting dark matter density distributions of dif- 
ferent Milky Way mass scale N-body simulations is at 
the level of a factor of two for distances larger than 
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^ 45 pc from the center of the galaxy as apparent in 
Fig. [T] For more information on dark matter from an 
astrophysical and pa rticle physics point of view see 
iBertone et aTl [2005j and references therein. 



2. The High Energy Stereoscopic System 

The High Energy Stereoscopic System (H.E.S.S.) 
is an array of Imaging Atmospheric Cherenkov Tele- 
scopes (lACTs) located in the Khomas Highland of 
Namibia at an altitude of ~ 1800 m above sea level. 
This study considers the H.E.S.S. I array of four iden- 
tical ~ 100 m^ mirror area lACTs operating in a 
square formation and optimized for the detection of 
~ 100 GeV to ~ 100 TeV 7-rays from astrophysi- 
cal sources. The instrumental field of view (FoV) of 
H.E.S.S. I is 5° in diameter. The energy and angu- 
lar resolution for the reconstruction of 7-rays is ^ 20% 
and 0.1° respectivly. For more technical info rma- 
tion on H.E.S.S. I see lAharonian et al. I |2008j | and 
references therein. 



2.1. Bacl<ground Subtraction and 
H.E.S.S. I Point Source Sensitivity 

The detection of 7-rays is complicated due to the 
necessity to distinguish 7-ray signal events from back- 
ground events from Cosmic Rays. H.E.S.S. I per- 
forms a discrimination between C osmic Ray and 7-ra y 



events in multiple steps at trigger ([Funk et al 
and subsequent data analysis ( Aharoni an et al 
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|2008j |) level. However, a perfect discrimination is not 
possible on an event by event basis and enforces the 
application of a background subtraction algorithm. 
Frequently used and well tested background subtrac- 
tion algorithm rely on the definition of a signal and a 
backg round region in the same FoV feee lBerge et al. I 
|2007j | for an overview) . The background region has to 
be constructed such that the instrumental acceptance 
in the signal and background region are equal. The 
number of expected background events in the signal 
region is then measured under the assumption that the 
Cosmic Ray TeV flux is isotropic. This background 
subtraction approach is well suited for the investiga- 
tion of signal regions which are much smaller than the 
FoV and leads typically to the ability to detect as- 
trophysical 7-ray point sources wit h a photon flux of 
~ 1% of the Crab nebula 7-ray flux ( Aharonian et al.~l 
[2008|) at a statistical significance of 5 standard de- 
viations within an observation time of ^ 25 h under 
optimal observation conditions for H.E.S.S. I. For 7- 
ray sources that are comparable in size or larger than 
the FoV of H.E.S.S. I the background region can how- 
ever not be constructed in this way and alternative 
methods have to be applied. 
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Figure 2: Signal (green) and background (red) region 
constructed with the rotated pixel background algorithm 
in galactic coordinates. Each pixel in the signal region is 
rotated out of the signal region at constant angular 
distance to the pointing direction. Known astrophysical 
7-ray sources and the galactic plane are excluded 
(yellow) . 



3. The Rotated Pixel Method 



The main technical problem in the search for a 7- 
ray signal from the galactic dark matter halo is that 
the halo extension is much larger than the H.E.S.S. 
I FoV. Traditional background subtraction algorithm 
that rely on the signal plus background and a back- 
ground measurement in the same FoV can thus not 
easily be applied. The H.E.S.S. collab oration has re- 
cently published (IH.E.S.S. Col. I \2Qll\ ] a search for a 
7-ray signal from self annihilating WIMPs in the cen- 
ter of the Milky Way dark matter halo which uses an 
adaptation of the r eflected regioii backg round subtrac- 
tion algorithm fsee lBerge et al. I |2007j | V The method 
employed relies on the steepness of the predicted dark 
matter density proflle in the vicinity of the galac- 
tic center and the assumption of the H.E.S.S. I in- 
strumental acceptance to be rotationally symmetric 
around the pointing position. A signal region of 1° 
angular radius around the galactic center is defined 
and subdivided into many pixel each of which is much 
smaller than the angular resolution of H.E.S.S. I. For a 
given observation run pointing towards a constant di- 
rection in celestial coordinates, each pixel in the over- 
lap of the signal region and the FoV is then rotated 
out of the signal region at constant angular distance 
to the pointing position to construct the background 
region for a run. Known astrophysical 7-ray sources, 
the galactic plane (|6| < 0.3°) and signal pixels for 
which no background pixel can be constructed for ge- 
ometric reasons are excluded from the analysis. Fig- 
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ure [T] shows the projected distances to the galactic 
center covered in the signal (magenta, tagged 'wobble 
on') and background (turquoise, tagged 'wobble off') 
region. The construction of the background region en- 
sures that the signal region is always closer in angular 
distance to the galactic center than the background 
region and the predicted average dark matter density 
in the signal region is thus larger than in the back- 
ground region. The background subtraction method 
relies thus not on a measurement of signal plus back- 
ground and background in the same FoV but on a 
measurement of signal plus background and 'less sig- 
nal' plus background in the same FoV. A search for 
a particle dark matter self annihilation signal is pos- 
sible by comparing the number of events detected in 
the signal region with the number of events detected 
in the background region. The H.E.S.S. collaboration 
did not detect a significant signal in a rotated pixel 
analysis of 112 h livetime of H.E.S.S. I observations of 
the galactic center region. The velocity averaged self 
annihilation cross section < av > was in consequence 
constrained to be smaller tha n ^ 3 • 10~^^ cm^/s f or 
WIMP masses around 1 TeV (|H.E.S.S. Col. I |201l| ). 
The derived limits on < cru > are currently the most 
constraining in the TeV WIMP mass range. 
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Figure 3: Simulated H.E.S.S. I events for the observation 
of a signal region with the fuU H.E.S.S. I FoV and two 
background regions in galactic coordinates. The pointing 
positions for the background regions have a symmetric 
offset in right ascension to the pointing position of the 
signal region. Known 7-ray sources and the galactic 
plane are excluded from the analysis and all exclusion 
regions are mutually applied to all observation regions in 
the FoV system. 



3.1. Limitations 

The application of the rotated pixel background 
subtraction algorithm limits the angular extension of 
the signal region by the demand to fit signal and back- 
ground region into one H.E.S.S. I FoV. A larger signal 
region could, depending on the signal to background 
ratio, increase the sensitivity of H.E.S.S. I to a WIMP 
self annihilation 7-ray flux from the galactic dark mat- 
ter halo. Without the constraint due to the rotated 
pixel background subtraction, the maximal signal re- 
gion could be increased to be complete H.E.S.S. I FoV. 
The maximal distance of the background region to 
the galactic center is also limited by the applied back- 
ground subtraction technique. If the background re- 
gion can be constructed much further away from the 
galactic center than possible with the rotated pixel 
technique, the expected 7-ray flux due to WIMP 
self annihilation in the background region would be 
smaller and the sensitivity would thus improve. 
In the following, two alternative background subtrac- 
tion techniques that rely on special observation strate- 
gies are described and it is investigated whether an in- 
crease in sensitivity to a WIMP self annihilation 7-ray 
flux from the Milky Way dark matter halo is possible. 



4. On/Off Bacltground Subtraction 

Figure [3] shows simulated H.E.S.S. I 7-ray events 
in three different FoVs. H.E.S.S. I extragalactic 'off 



data', i.e. observation runs with no 7-ray source in 
the FoV, performed at a zenith angle of ~ 20° were 
used to obtain a realistic 7-ray candidate event rate af- 
ter H.E.S.S. I standard 7-ray selection criteria (cuts). 
Additionally, the dependence of the 7-ray candidate 
event rate after cuts as a function of the offset to the 
observation position ('radial acceptance') and the en- 
ergy dependence of the event rate after standard cuts 
together with the instrumental dead time as obtained 
from extragalactic off data are used to simulate real- 
istic observations. Figure |3] shows in the center the 
simulated signal region events in the vicinity of the 
galactic center distributed over a complete H.E.S.S. 
I FoV. The observation of two background regions 
whose observation positions have a symmetric offset of 
±35 min in right ascension to the signal region obser- 
vation position is simulated. Known 7-ray sources and 
the galactic plane (|&| < 0.3°) are excluded from the 
analysis and all exclusion regions are mutually applied 
to all observed FoVs respectively. For each simulated 
FoV, 50 h observation livetime were generated. 
The simulation models a peculiar observation strat- 
egy where always three H.E.S.S. I standard observa- 
tion runs are taken consecutively such that the zenith 
and azimuth pointing angle range is always equal for 
a sequence of three runs. One observation sequence 
starts with the observation of a background region 
with —35 min offset in right ascension from the sig- 
nal region pointing position and is performed with a 
standard length of 33 min. The array pointing moves 
forward by 35 min in right ascension immediately af- 
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ter the end of the first observation within a realistic 
run transition time of 2 min. The second run of the 
sequence observes then the signal region for 33 min. 
The third observation is again scheduled for 33 min 
immediately after the second observation by moving 
the array pointing again by 35 min forward in right 
ascension within 2 min. Finally, the number of events 
after the application of standard 7-ray cuts in the 
signal and and the corresponding number of events 
in background region are compared. The observation 
strategy ensures as stated that the zenith and azimuth 
pointing angle range is always the same for every run 
within a run sequence. The strong dependence of the 
H.E.S.S. I acceptance on the pointing zenith and az- 
imuth angle does in consequence not influence the cal- 
culation of a 7-ray excess as the dependence is the 
same for the signal and the background region. Re- 
peated taking of three run sequences in this obser- 
vation pattern can accumulate the simulated amount 
of data. The advantage is that obviously the signal 
region can be as large as the H.E.S.S. FoV (see also 
the blue region in fig. [T]) and the background regions 
(see also the yellow region in fig. [T]) are further away 
from the signal region with the result of less expected 
WIMP self annihilation 7-ray flux in the background 
region. This leads to the expectation that the sen- 
sitivity of the On/Off observations to a WIMP self 
annihilation 7-ray flux in the Milky Way halo is bet- 
ter than for the rotated pixel background subtraction 
technique. On the other hand, the method needs 2/3 
of the total observation time for the estimation of the 
expected number of background events in the signal 
region where the rotated pixel algorithm can use the 
complete observation time. It is a non trivial ques- 
tion which of the two competing factors is eventually 
dominant and therefore this realistic simulation is per- 
formed. 



5. Driftscan Observations 

Figure |4] shows the simulation of events recorded 
in driftscan observations of the galactic center region. 
Driftscan runs are scheduled at a constant zenith and 
azimuth angle array pointing. The celestial pointing 
at the beginning of a run is chosen such that the galac- 
tic center region 'drifts through' the FoV after approx- 
imately half of the 68 min observation time. As for 
the simulation of On/Off events, realistic event rates 
for the observation at ^ 20° zenith angle after stan- 
dard H.E.S.S. I 7-ray cuts have been used. Also the 
radial acceptance as well as the energy dependence of 
the events after 7-ray cuts and the instrumental dead 
time resemble real observation conditions. A total of 
150 h of observation time is simulated. The observed 
region is divided into many pixels each of which is 
much smaller than the H.E.S.S. I angular resolution. 
Each run is divided in time into one signal region en- 
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Figure 4: Simulated H.E.S.S. I events for the driftscan 
observation of the Milky Way dark matter halo in 
galactic coordinates. Simulated is the repeated taking of 
data at constant zenith and azimuth pointing where the 
galactic center region 'drifts' through the FoV after 
approximately half of the observation time. Each run is 
subsequently divided into two parts, the signal region 
enclosing the galactic center region in the lower panel 
and the background region in the upper panel. Known 
7-ray sources and the galactic plane (|&| < 0.3°) are 
excluded. Exlcuded regions are treated with a special 
method (see text) to guarantee the same instrumental 
acceptance in the signal and background region. 



closing the galactic center and one background region. 
The division of a driftscan run into two halves is re- 
alized such that the product of the average time a 
pixel stays in the H.E.S.S. I FoV with the total solid 
angle of the region is the same for the signal and back- 
ground region. As for the rotated pixel technique and 
the On/Off observations all known 7-ray sources and 
the galactic plane (|6| < 0.3°) are excluded from the 
analysis. Excluded pixel are shifted mutually in right 
ascension between signal and background region at 
constant declination and special care is taken to guar- 
antee that the shifting of the exclusion regions does 
not imbalance the instrumental acceptance of signal 
and background region. Every pixel excluded in the 
signal region is shifted along right ascension and con- 
stant declination as close as possible to the edge of 
the background region to minimize the expected 7- 
ray flux from WIMP annihilation in the background 
region. Similarly, every pixel excluded in the back- 
ground region is shifted as close as possible to the 
edge of the signal region to maximize the expected 7- 
ray flux from WIMP annihilation in the signal region. 
The observation strategy and the subsequent data 
analysis ensure again, as for the On/Off and rotated 
pixel method, an equal instrumental acceptance in the 
signal and background region. Also similar to the 
On/Off method, the signal region is larger than for 
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Figure 5: Sensitivity comparison between the rotated 
pixel method (blue), the On/Off method (green) and the 
driftscan method (red). Shown is the sensitivity (average 
95% C.L. upper limit) to the velocity averaged self 
annihilation cross section < cru > as a function of the 
assumed WIMP mass in case that no signal can be 
detected. 



the rotated pixel method and the background region 
is further away from the galactic center than possible 
for the rotated pixel method. These effects should in- 
crease the sensitivity to a WIMP annihilation 7-ray 
flux from the galactic halo. On the other hand, the 
increase in sensitivity competes with the decrease in 
sensitivity due to observation time only used for back- 
ground observations and the average expected signal 
flux per signal pixel is smaller than for the On/Off 
method because the signal FoV is larger. Again, it is 
non trivial to estimate which effect dominates and a 
realistic simulation that compares all three discussed 
methods under equal conditions is necessary. 



6. Sensitivity Comparison 

Figure [5] shows the sensitivity of the rotated pixel, 
the On/Off and the driftscan method to the veloc- 
ity averaged self annihilation cross section < av > 
as a function of the WIMP mass. The Einasto 
parametrization for the dark matter density distri- 
bution of t he Milky Way is a ssumed, parameters are 
taken from lPieri et al. 1 [2011 1. The continuous 7-ray 
spectrum exp ected from the WIMP self ann ihilation is 
adopted from lTasitsiomi and Olinto ] |2002l |. For each 
investigated method, the On/Off, the driftscan and 
the rotated pixel method a total of 150 h of observa- 
tion time is simulated. No WIMP annihilation 7-ray 
signal is simulated but only background events such 
that no significant 7-ray excess is resulting in all cases 
and an upper limit on < av > can be derived which 
serves as a sensitivity estimation for the method un- 



der consideration. The upper limit on < crw > itself 
is calculated at 95% confidence level using H.E.S.S. I 
effective areas a nd energy threshold look ups for stan- 
dard Hillas cuts ( Aharonian et al.~l (2008j |) and with a 
method that leads to res ults that are com patible with 
the method employed in IH.E.S.S Col. I 2011] under 
equal conditions. 

Figure [5] shows that the sensitivity to < crw > of 
each method investigated is comparable in the order 
of magnitude. The On/Off method is slightly more 
sensitive (~ 20% at 1 TeV WIMP mass) than the ro- 
tated pixel method. The driftscan method is less sen- 
sitive a factor of 2 at 1 TeV WIMP mass) than the 
rotated pixel method and the On/Off method. How- 
ever, the sensitivity of the driftscan method can be 
increased compared to the result shown in Fig. [5] 
when the size of the driftscan background region is 
increased at the cost of the size of the signal region. 
In an optimized signal to background region size case 
the sensitivity of the driftscan method agrees with the 
sensitivity of the On/Off and rotated pixel method 
within - 30%. 

Overall it can be stated that the increase in sensi- 
tivity of background subtraction methods alternative 
to the rotated pixel method is not as large as could 
be expected. The reason for this is the large amount 
of observation time that the alternative background 
subtraction methods need only for background obser- 
vations. On the other hand, it is obvious that the al- 
ternative background subtraction methods can be sen- 
sitive to a WIMP annihilation 7-ray flux even in the 
case when the dark matter density profile is for some 
reason nearly constant in the vicinity of the Milky 
Way halo center but only decreasing at large distances 
to the galactic center. Additionally, the sensitivity of 
the driftscan and On/ Off method is larger than the 
sensitivity of the rotated pixel method if the center of 
the Milky Way dark matter halo is not spatially coin- 
cident with the galactic center but off set by 1° — 2° 
as possibly indicated in iKuhlen et al I (2012j . In con- 
clusion it is (despite of the sensitivity to the discussed 
benchmark model with an Einasto dark matter den- 
sity profile being comparable to the rotated pixel tech- 
nique) of interest to record at least a limited amount 
of observational data with the discussed observation 
strategies. 
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